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Abstract: The status of the pair-production dip as a spectral feature, produced by interaction 
of Ultra High Energy extragalactic protons with CMB is discussed. 



Greisen-Zatsepin-Kuzmin (GZK) cutoff pQ is 
the most spectacular prediction for Ultra 
High Energy Cosmic Ray (UHECR) spectrum, 
which status is still uncertain in the present ob- 
servations. As physics is concerned, detection 
of the GZK cutoff means discovery of UHE pro- 
ton interaction with CMB radiation. Another 
prediction for interaction of UHE protons with 
CMB is pair-production dip, the spectral fea- 
ture originated from electron-positron pair pro- 
duction by extragalactic UHE protons propa- 
gating through CMB: p + 7cMB ^ p + e + + e~. 
Originally proposed for diffuse spectrum in 
early work [2], this feature has been studied 
recently in Refs. [H [4]. An alternative ex- 
planation of the observed pair-production dip, 
widely discussed now [5], was first put for- 
ward in works 6 a and [7] in terms of a two- 
component model as the transition from galac- 
tic to extragalactic cosmic rays. Being a quite 
faint feature, the e + e _ -production dip is not 
seen well in the naturally presented spectrum 
log J(E) vs. log E. The dip is more pronounced 
when analyzed in terms of the modification 
factor [111], T)(E) = J p (E)/J£ nm (E), where 
J P (E) is the spectrum calculated with all en- 
ergy losses included, and J™ m (i?) is the un- 
modified spectrum calculated with adiabatic 
energy losses only. The observed modifica- 
tion factor is given by T^bs oc J y, s (E) / E~ la , 
where J bs(E) is the observed spectrum and 
7 g is the exponent of the generation spectrum 



Qgon(Eg) oc Eg 7s in terms of initial proton en- 
ergies Eg. 

The pair-production dip is clearly seen in the 
energy-dependence of n(E) and is reliably con- 
firmed O IU[9] by observational data, as Fig.Q] 
shows. The comparison of the predicted dip 
with observational data includes only two free 
parameters: exponent of the power-law gener- 
ation spectrum j g (the best fit corresponds to 
7 g = 2.6 — 2.7) and normalization constant to 
fit the e + e~-production dip to the measured 
flux. The number of energy bins in the differ- 
ent experiments is 20 - 22. The fit is charac- 
terized by x 2 /d.o.f. = 1.0 - 1.2 for AGASA, 
HiRes and Yakutsk data. For the Auger data 
% 2 /d.o.f. is larger mainly due to the low flux in 
the first energy bin at E w 45 EeV where mea- 
surements are made with the help of surface 
detectors (see Fig. [T]). 

The theoretical pair-production dip has two 
flattenings: one at energy Ej, » 1 X 10 18 eV 
and the other at E a « 1 x 10 19 eV. One can see 
that &t E < E}, the experimental modification 
factor, as measured by Akeno and HiRes, ex- 
ceeds the theoretical modification factor. Since 
by definition modification factor must be less 
than one, this excess signals the appearance of 
a new component of cosmic rays at E < Ef, = 
1 x 10 18 eV, and thus the transition from ex- 
tragalactic to galactic cosmic rays, starting at 
energy E b . 

The second flattening automatically explains 
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Figure 1: The predicted pair-production dip in comparison with Akeno-AGASA, HiRes, Yakutsk 
and Auger data [10]. The first three experiments confirm dip with good % 2 /d.o.f. « 1.0 — 1.2, while 
the Auger data are characterized by larger x 2 /d.o.f., mainly due to the bin at energy E as 45 EeV. 
The data of Fly's Eye [10] (not presented here) confirm the dip as AGASA, HiRes and Yakutsk 
detectors do. 



the ankle, the feature seen in all experiments 
starting from Haverah Park in the end of 70s. 

The position and shape of the dip is robustly 
fixed by interaction with CMB and can be used 
for energy calibration of the detectors. 

The systematic errors in energy measurements 
are high, from 15% in AGASA to 22% in 
Auger. To calibrate each detector we shift the 
energies by factor A to reach minimum y 2 in 
comparison with theoretical dip. We obtain 
these factors as A^ = 0.9, Ay Q = 0.75 and 
\ m = 1-2 for AGASA, Yakutsk and HiRes de- 
tectors, respectively. Recently, AGASA collab- 
oration has reduced their energies by 10% in- 
deed, based on reconsideration of energy deter- 
mination. After energy calibration the fluxes 
given by AGASA, HiRes and Yakutsk detec- 
tors agree with each other in most precise way 



(see Fig. [2]). The Auger flux is noticeably be- 
low the flux shown in Fig. O 

Concerning the calibration two remarks are in 
order. 

i) After calibration the discrepancy between 
AGASA and HiRes data at the highest ener- 
gies diminishes to the level of 2.5 a, but the 
AGASA excess over the theoretical flux with 
the GZK cutoff remains statistically signifi- 
cant. The better agreement between highest 
energy AGASA and HiRes data implies some 
trial theoretical spectrum between AGASA 
and HiRes data. 

ii) One can see that calibration with help of the 
pair-production dip implies decreasing energies 
measured by on-ground methods (Xa = 0.9 
and Xy a = 0.75) and increasing the energies 
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Figure 2: The fluxes from Akeno-AGASA, HiRes 
by the e + e~ -production dip. 

measured by fluorescent method (Xm — 1-2). 
It might be considered as an indication to the 
difference in measuring energies by these two 
methods. The predicted shape of the e + e~- 
production dip is quite robust [H [9]: it is 
modified very weakly when the discreteness 
in the source distribution and their inhomo- 
geneities are taken into account, and differ- 
ent regimes of propagation (from rectilinear 
to diffusive) are considered. The cosmologi- 
cal evolution of the sources, e.g. with parame- 
ters inspired by observations of active galac- 
tic nuclei, also results in the same shape of 
the dip. The pair-production dip is modified 
strongly when the fraction of nuclei heavier 
than protons is high at injection, which im- 
poses some restrictions to the mechanisms of 
acceleration operating in UHECR sources [5]. 
The shape of acceleration spectrum needed for 
the e + e~ -production dip agrees with standard 
ones 7 S = 2 for non-relativistic shock accelera- 
tion or 7 g = 2.2 — 2.3 for relativistic shock. The 
effective 7 9 = 2.6 — 2.7 needed at ultra high en- 
ergy is naturally explained by distribution of 
sources over maximum energy of acceleration 
or luminosity [USE]. 

On the basis of the predicted dip and the cal- 
ibrated data of AGASA, HiRes and Yakutsk 
detectors we can plot the spectrum and flux 
in the energy region 1 x 10 18 — 1 x 10 20 eV 
as it is shown in Fig. [3J In the energy interval 
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Yakutsk detectors before and after calibration 



(0.1 — 10) x 10 19 eV the theoretical uncertainties 
in the predicted spectrum are relatively small 
and are mainly given by uncertainties in dis- 
tances between sources. These uncertainties 
dramatically increase at E > 1 X 10 20 eV. 

In Fig. [3J the spectra are shown for proton- 
dominated flux with distances between sources 
in the range (1 — 60) Mpc. Therefore the be- 
ginning of the GZK cutoff in the energy range 
(5 — 10) x 10 19 eV is predicted in the dip-based 
model with small uncertainties. At larger en- 
ergies the spectrum of GZK feature is very 
model dependent: apart from distances be- 
tween sources it depends on fluctuations in lu- 
minosities of the nearby sources, in distances 
between them, and by maximum acceleration 
energy E max (see [3J for calculations). One can 
see from Figs.Q]and[3jthat the beginning of the 
GZK cutoff in energy range (5 — 10) x 10 19 eV is 
confirmed by all detectors, including AGASA. 
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Figure 3: The dip-based prediction for diffuse spectrum in energy range f x fO 18 — fx fO 20 eV. 
The calculated dip is normalized by calibrated AGASA- Yakutsk data as shown in Fig. [2] The 
calculated diffuse energy spectrum is displayed for different distances d between sources in the range 
f — 60 Mpc. This presents the largest theoretical uncertainties in energy range (1 — 8) x 10 19 eV. The 
both uncertainties in spectrum in the interval (1 — 8) x 10 19 eV due to the dip-based calculations 
and measurements by AGASA, Hires and Yakutsk detectors are small enough, and the beginning 
of the GZK cutoff at (5 - f 0) x f 19 eV is reliably predicted. 
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